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ABSTRACT 

We study free string propagation in families of plane wave geometries developing strong 
scale-invariant singularities in certain limits. We relate the singular limit of the evolution for 
all excited string modes to that of the center-of-mass motion (the latter existing for discrete 
values of the overall plane wave profile normalization). Requiring that the entire excitation 
energy of the string should be finite turns out to be quite restrictive and essentially excludes 
consistent propagation across the singularity, unless dimensionful scales are introduced at 
the singular locus (in an otherwise scale-invariant space-time). 
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1 Introduction 



String propagation in strong gravitational waves has attracted a considerable amount of 
attention on account of a few highly special properties of such space-times (see [H E] among 
other publications). For one thing, the structure of the curvature tensor in plane gravita- 
tional waves implies that these solutions to Einstein's equations (coupled to appropriate 
matter fields, if necessary) remain uncorrected lU [3] in a number of higher derivative exten- 
sions of general relativity (and, in particular, they do not receive any a'-corrections when 
introduced as backgrounds in perturbative string theories). Furthermore, the correspond- 
ing light cone Hamiltonian of string a-models turns out to be quadratic and admits a fairly 
thorough analytic treatment. (This class of backgrounds also admits a natural formulation 
of the matrix theory description of quantum gravity [HE].) 

In this publication, we shall concentrate on backgrounds of the following form (so-called 
exact plane waves): 

d d 

ds^ = -2dX+dX~ - F(X+) ^(X')2(rfX+)2 + Y^idX' f- (1) 

i=l i=l 

This representation of the metric is often called the Brinkmann form. The case of constant 
F(X+) corresponds to supersymmetric plane waves studied in [6], and it is quite different 
from the the rapidly varying F{X~^) we intend to consider. (A coordinate transformation 
can be performed into the so-called Rosen coordinates eliminating the dependence of the 
metric on the transverse coordinates X*. The resulting metric depends on X+ only and 
displays manifestly a plane-fronted space-time wave propagating at the speed of light. How- 
ever, the Rosen parametrization tends to suffer from coordinate singularities, and we shall 
work with the Brinkmann form.) 

The function F(X+) contained in ([T]) is completely arbitrary, and one may ask, for 
example, what happens to quantum strings propagating in such space-times when the wave 
profile F{X~^) develops an isolated singularity. This question is of some interest per se, 
since studies of string theory in the presence of space-time singularities have played a 
pivotal role in the development of the subject (and, in this particular case, we are dealing 
with singularities in time-dependent backgrounds). Additional heuristic justification for 
our studies is provided by the observation that plane waves of the type ([T]) with 

— (2) 

arise as Penrose limits of a broad class [7J of space-time singularities (including the Friedmann- 
Lemaitre- Robertson- Walker cosmological singularities). With F{X~^) of ([2]), the metric ([T]) 
is invariant under scaling transformations X"*" aX"*", X^ X^ /a (identical to Lorentz 
boosts in flat space-time). Note that this type of singularities is considerably stronger than 
the so-called "weak singularities" of [8]. 

Free string propagation on ([1]) with F{X~^) given by ([2]) has been previously studied 
in [2]. In particular, it was suggested in that publication that the question of propaga- 
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tion across the singularity in the metric can be addressed by employing analytic 

continuation in the complex X^-plane. We believe that this issue merits further elucidation. 

In the context of string theory and related approaches to quantum gravity, there is a 
general expectation that the space-time background used for formulating the theory should 
satisfy some stringent consistency conditions. For perturbative string theories, these con- 
ditions take the form of the appropriate supergravity equations of motion together with an 
infinite tower of a'-corrections. For non-singular plane waves, all the a'-corrections vanish 
automatically on account of the special properties of the Riemann tensor corresponding to 
these space-times. For singular space-times, the question of background consistency condi- 
tions at the singular point appears to be extremely subtle. Indeed, what should replace the 
supergravity equations of motion at the singular point where they obviously break down? 
Ad hoc prescriptions are not likely to produce meaningful results under these circumstances. 

One approach to formulating string theory in backgrounds ([T]|2]) is to resolve the singular 
plane wave profile into a non-singular function, perform the necessary computations and 
see if the result has a meaningful singular limit. (This approach was advocated in [IJ, where 
a conjecture was made that for certain choices of the plane wave profile, taking a singular 
limit may result in well-defined transition amplitudes. We intend to consider this question 
quantitatively.) Note that, for the resolved space-times of this sort, perturbative string 
background consistency conditions are automatically satisfied to all orders in a'. The only 
non-trivial question is the existence of a singular limit. 

But how should one resolve? We want to construct a function F(X+, e) in such a way 
that 



everywhere away from X'^ = 0. There is in principle a large amount of ambiguity associated 
with such resolutions. One class appears to be very special however. The background ([T]l2]) 
possesses a scaling symmetry and does not depend on any dimensionful parameters. It is 
natural to demand that this symmetry should be recovered when the resolution is removed. 
This will happen if the resolved profile F{X~^, e) does not depend on any dimensionful 
parameters other than the resolution parameter e. In this case, on dimensional grounds. 



for large values of rj, with some b > 2. Note that the fact that the original background 
possesses a certain symmetry (away from = 0!) in no way implies that we must resolve 
in a way consistent with this symmetry. For resolved profiles different from (jlj), the limit 
of the metric may still be given by ([3]) away from X~^ = (and thus be scale invariant), 
but additional dimensionful scales may become buried inside the singularity at X~^ = 
(in a way that only affects processes involving singularity crossing). One would need some 




(3) 




(4) 



The limit ([3]) will be recovered if 




(5) 
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strong physical rationale for introducing such scales buried at the singular locus, and in the 
present publication we shall simply study the "scale- invariant" resolutions (jlj). 

The structure of the paper is as follows: we will first derive the Hamiltonian for a 
free string in the background ([I112D- Then we recapitulate the main results of [2] for the 
evolution of the center-of-mass motion across the plane wave singularity. We extend this 
analysis to the evolution of excited string modes. We conclude by discussing stringent 
conditions arising if one demands the total mass of the string to remain finite after it 
crosses the singularity. 



2 Free strings in plane waves 

Due to the presence of covariantly constant null vectors in plane wave geometries, the 
string theory cr-model can be analyzed explicitly in such backgrounds, and reduces to a set 
of independent classical time-dependent harmonic oscillators. In this section, we re-state 
this familiar material in a way convenient for our present investigations. 

2.1 The light cone gauge 

String worldsheet fermions are free in plane wave backgrounds [lOj. We shall therefore 
concentrate on the bosonic part of the string action, given by 



(6) 



47ra 

We choose light-cone gauge = a'p^r and gauge-fix the metric, 

dei{gab) = -1, d^gaa = 0, (7) 
to obtain the following Lagrangian, where we have solved for Qrr'- 



'^'^^^\{a'p+r/t){X'f 
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2gra [a'p+d^X- - drX^d^X^) + g-l{l - gl,) Y,{d^X^f - \a'R^'^^ . 



We rescale e = e'a'p"^. 



= -^^{r/e'l (9) 



g2 y I J 



and from here on, we will denote worldsheet time t = t and write e instead of e'. The 
(T-dependent part of the oscillator X~ is non-dynamical and enforces gra = 0. The cr- 
independent part of the oscillator X~ can be eliminated as a constraint {gcra = 1), the 
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(dynamically non-trivial) coupling to the dilaton disappears (see, e.g., [2]), and we can 
write the following worldsheet Hamiltonian 



i=l ^ ^ 



(10) 



where Pi are momenta conjugate to X*. We will now choose units in which a' = 1. If we 
Fourier transform the a-coordinate, 

X\t, a) = X^(t) + (cos (na) X;(t) + sin (na) X^t)') , (11) 

n>0 

we obtain a set of time-dependent harmonic oscillator Hamiltonians 

oo d 

H = Y,Y.Hm, (12) 
n=0 i=l 

^^oi = ^ + 4n(*A)^. (is) 



. (^y!±(M , („. , W)l±™!. ,14) 

2.2 WKB solution for time-dependent harmonic oscillator 



The Hamiltonian (11211 is quadratic and the solution to the corresponding Schrodinger equa- 
tion, 



(15) 



can be found using WKB techniques, which are exact for quadratic Hamiltonians. From 
(fTSD it follows that 

if we separate variables as 

8 

0(t;x)=nn<^"(^'^-)- (1^) 

n i=l 

We then take the WKB ansatz 

0^(t;X) = A(ti,t)exp (25'eZ;„[X{_„,ti|X;,t]) , (18) 

where S'c/;„[XJ „, ti|X^, t] is the "classical action" evaluated for the path going from X{ ,j at 
the time ti to at the time t, 

s4xir.Mx:,t] = j\t' -{f + T2^i'-\]^^\- (19) 



li An{ti,t) satisfies 

-2^^An{h,t) = Mtut)^^S,i[Xi^,t,\Xlt], (20) 

then (fTSll satisfies the original Schrodinger equation exactly. 

Up to normalization, a basis of solutions, labelled by the initial condition X^(ti) = Xl^, 
is given by [9], 

^^^-^ ^n) ~ n |j [ W,n)'5t.C - (X;)29,,C + 2X1^X1]^ , (21) 

where C(ti,t2) (suppressing the index n) is a solution to the "classical equation of motion" 
for the time- dependent harmonic oscillator Hamiltonian 

dlC{hM) + [n" + ^^^(tsA)^ C{hM) = 0, (22) 

with initial conditions specified as 

C{h,h)\t,=t, = 0, dt,C{hM)\t,=t, = 1. (23) 

We shall refer to C{ti,t2) as "compression factor", since it describes convergence of solutions 
to the corresponding harmonic oscillator equation starting at the same point at the moment 
ti- (If C(ti, ^2) vanishes, then t2 is a focal point, as the difference between any two solutions 
with the same initial position X{ti) is proportional to C{ti,t2)-) A useful representation of 
C(ti, is given by 

... f{h)h{t2) - f{t2)h{t,) 

C[hM)- , (24) 

where f{t) and h{t) are two independent solutions to the differential equation under con- 
sideration, and the Wronskian W is given by 

W[f,h] = fh-hf. (25) 

To derive the singular limit of the wavefunction ( l2Ti) it is sufficient to study the singular 
hmit of (EMI. 



3 The singular limit for the center-of-mass motion 

For the n = mode, we obtain as the "classical equation of motion" 

X + ^^{t/e)X = 0. (26) 

We need to study the e ^ limit of the solution that obeys the initial conditions 

X(ti) = 0, X(ti) = l, ti<0. (27) 
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The singular limit of solutions to this equation has been analyzed in [9]. Performing 
a scale transformation Y{r]) = X{rie), with t] = t/e, removes the e-dependence from the 
equation, leaving 

—Y + \n{ri)Y = 0. (28) 

This scale transformation is possible because our initial singular metric was scale-invariant 
and we have resolved it as in (jlj) without introducing any dimensionful parameters besides 
e. The existence of a singular limit is then translated [9] into constraints on on the asymp- 
totic behavior of solutions to (l28l) . These "boundary conditions at infinity" are strongly 
reminiscent of a Sturm-Liouville problem, and it is natural that a discrete spectrum of A is 
singled out by imposing the existence of a singular limit. 

For the specific asymptotics of our resolved profile ([5]), it can be shown |9] that, in the 
infinite past and infinite future, the solutions approach a linear combination of two powers 
(denoted below a and 1 — a, with a being a function of kX, cf. dHIH])). This power law 
behavior simply corresponds to the regime when the second term on the right hand side 
of ([5]) can be neglected compared to the first. It is then convenient to form two bases of 
solutions, one asymptotically approaching the two powers (dominant and subdominant) at 

7] -OO, 



Yi-iv) = \vr+o{\vn, 

and another behaving similarly at ^ +oo 

Y,4v) = \vr+oi\vr), 



Y,4r^) = Ir^l'-'^- + o{\r^\'- 



Y2+{v) = \v\'-''^ + o{\v\'-^+), 



(29) 



(30) 



where a± is given by 



a± 



Xk 



(31) 



(We are temporarily assuming that k can take two different values k± for the positive and 
negative time asymptotics, a possibility that will be discarded shortly.) The two bases need, 
of course, to be related by a linear transformation: 



Y2-{V) 



(32) 



where Q{X) is a 2 x 2 matrix whose determinant is constrained by Wronskian conservation 
as 

W[Yi^, Y2-] = W[Yi+, Y2+] det Q. (33) 

The singular limit has been rigorously considered in |9j, but the results can be un- 
derstood heuristically from the following argument. Imagine one is trying to construct a 
solution Y to (l28l) satisfying some (e-independent) initial conditions at ?7i = ti/e < 0. This 
solution can be expressed in terms of Yi_ and (a complete basis) as 



Y = CiFi_ + C2F2- 



(34) 



6 



Since the initial conditions are specified at rji = ti/e, the asymptotic expansions fl29l) are 
vahd. There needs to be a non-trivial contribution from both Yi_ and Y2- in the above 
formula in order to satisfy general initial conditions. Hence, the two terms on the right 
hand side should be of order 1. Therefore, we should have 

Ci=0(e»-), C2 = 0(ei-»-). (35) 

If we now apply fl^ and fl5U]) to evaluate at a large positive rj = t2/e, the powers of e in 
Ci and C2 will combine with the powers of e originating from Yi+ and ^2+ and yield 

Yih/e) = gn(A)CO(e'^-«+) + Qi2(A)4-"+0(e'^-+'^+-i) ^ ^ 

(36) 

+Q2iix)tToie'-"^) + g22(A)4"'^^o(6"+-'^-). 

Since a+ and a_ are greater than 1/2, this expression can only have an e — limit if 
a-|_ = a_ (i.e., fc+ = A;_ and we can set both equal to 1 by redefining A) and (521(A) = 0. 
The latter condition implies that the absolute normalization A of the plane wave profile fl{rj) 
will generically lie in a discrete spectrum, dependent on the specific way the singularity is 
resolved, i.e., the shape of A particular exactly solvable example for this discrete 

spectrum (there called "light-like reflector plane") has been given in [Hj. With Q2i(A) = 
and det Q = —1, the matrix Q can be written as 



Q 



q q 
-1/q 



(37) 



with q being a real nonzero number {q does not affect the singular limit). For fiat space- 
time we have q = 1 and for the "light-like reflector plane" of [TT] we have q = —1. In the 
singular limit, a basis of solutions is given by 



Y,{t) = i-tr, Y^it) = (-t)i-», t < 0, 

1 

-i 
Q 



Yi{t)=qt\ Y2{t) = --t^-\ t>0. (3^ 



4 The singular limit for excited string modes 

Following our general discussion of free strings in plane wave backgrounds, the evolution of 
excited string modes is described by time-dependent harmonic oscillator equations 

^X(t) + [n" + j,n{tle)^ X{t) = 0. (39) 

Solutions for the wavefunctions of the excited string modes can be expressed in terms of 
a particular solution to this equation C(ti,t2) defined by (12211231) . Hence, to analyze the 
singular (e —>■ 0) limit of the excited modes dynamics, it should suffice to analyze the 
singular limit of C{ti,t2). Because is finite, it is natural to expect that it does not affect 
the existence of the singular limit (governed by the singularity emerging from Q{t/e)). We 
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shall prove that it is indeed the case for positive A (for negative A unstable motion of the 
inverted harmonic oscillator leads to divergence^. 

To derive C(ti,t2) for equation (|39|) we use the following strategy: the differential equa- 
tion fl5^ is linear and any solution X{t2) ai t = t2 can be written in terms of a "transfer 
matrix" T that only depends on the initial and final times, 



The transfer matrix can be expressed as 



T{ti,t2 



(40) 



-dtfi{ti,t2) C{ti,t2) 



(41) 



where c^f. and dt^ indicate differentiation with respect to the first and second argument 
respectively. The transfer matrix is completely determined once C(ti,t2) has been deter- 
mined, and vice versa. We will now use the fact that transfer matrices of subintervals 
are combined by ordinary matrix multiplication. Dividing the solution region into three 
sub-intervals, we shall calculate the transfer matrices for each sub-interval k and apply 
multiplication to construct the total transfer matrix. The sub-intervals shall be chosen as 
indicated in the following figure: 

I II III 







t2 



We use to indicate a time that will approach zero in the singular limit as 

l—c7c 



(42) 



with t staying finite in relation to the "moments of observation" ti and t2- The number 
c (between and 1) will be chosen later as needed for our proof. On each interval, we 
can write the transfer matrix in terms of the "compression factor" C^. Using matrix 
multiplication to construct the full transfer matrix T, we can now deduce an expression for 
the "compression factor" of the complete interval 

C(ti, t2) =Ciih, -QdtAii-te, QduCiiiiU, t2) - dt,Ci{ti, ~QCn{-te, U)dtfiin{U, ta) 

(43) 

in terms of the "compression factors" of the three sub-intervals. Once again, (9^. and dt^ 
differentiate C with respect to its first and second argument (initial and final time). 



^More specifically, the divergences arise from sub-leading infinities in the position of the inverted har- 
monic oscillator, while the leading infinities cancel. (Such sub-leading infinities are absent for the center- 
of-mass motion analyzed in section [3l hence no analogous divergences in that case.) Further details are 
given in section 4.4 and appendix A. 
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To study the existence of the singular hmit of C{ti,t2), we will use the following strat- 
egy: for two linear differential equations related by a small perturbation we will establish 
a bound on the difference between perturbed and unperturbed solutions with the same 
initial conditions. This bound will, of course, apply to Ck- For each of the three sub- 
intervals introduced above, we will consider a simplified differential equation that is a good 
approximation to equation (l39l) on the corresponding interval: 

• Region I and III: X{t) + (n^ + A/t^) X(t) = (related to Bessel's equation); 

• Region II: X{t) + X/e'^Q{t/e)X{t) = (equation of motion for the zero mode). 

Then, on each sub-interval, Ck can be written as the sum of a simplified "compression 
factor" Ck satisfying the simplified differential equation on this sub-interval, plus a small 
perturbation 6Ck- We will prove that, in the singular limit, the 6Ck will drop out of the 
expression for the total "compression factor" C{ti,t2). 

Most of this section is dedicated to implementing the proof we have just outlined. The 
reader primarily interested in the discussion of the singular limit and content with the 
general sketch given above can skip to section 14. 4[ 

4.1 Bounds on solutions to perturbed differential equations 

In view of the subsequent application to the singular limit analysis, we would like to bound 
the difference 6X between the solution X{t) of a perturbed differential equation, 

— X + (T + 5T)X = 0, (44) 

and the solution X{t) of an unperturbed differential equation, 

§^X + TX = 0, (45) 

where we take 

X = X + SX (46) 
and demand that the initial conditions remain unchanged: 

X{to) = X(to), dtX{to) = dtX{to). (47) 



If we substitute (1461) and (1451) into ( l44l) we obtain a differential equation for the perturbation 
on the solution. ^ 

^6X + T{t)6X = -6T{t) (X + 6X) . (48) 
A formal solution to (H8l) is given by 



/oo 
Gr{t, t')ST{t') (X(t') + SX{t')) dt', (49) 
-oo 
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with the Green function Gr{t,t') satisfying 



+ T(t)^ Gr{t,t')=5{t~t') 



(50) 



and initial conditions 

Gr{t,t%=to = 0, dtGrM\t=t, = 0. (51) 

Therefore, we can write the Green function in terms of the "compression factor" C of the 
unperturbed equation fH5l) . where C obeys the same initial conditions as in fl23|) : 

G.(M') (52) 
lO otherwise. 

To obtain a bound on 5X we will invoke the so-called Gronwall inequality 



4.1.1 The Gronwall inequality 

Let I = [A,B]. Assume (3 and a real valued and continuous on I and /9 > 0. If u is 
continuous, real valued on I and satisfies the integral inequality 

u{t) < a{t) + f f3{s)u{s)ds, t e I, (53) 

J A 

then ^ ^ 

u{t) < a{t) + j (3{s)a{s) exp (^j ^ (3{r)dr^ ds, t G /. (54) 

Proof: First we define ^ 

z{t) = / f3{s)u{s)ds, t e I. (55) 

J A 

Then, after differentiation and using the initial assumption fl53|) . we obtain 

z'{t) = (3{t)u{t) < (3{t)a{t) + f3{t)z{t). (56) 
Using the line above we write 



exp J (3{u)duj z(s)J = exp j^- J (3{r)dr j {z'{s) - l3{s)z{s)) (57) 

< (3{s)a{s)exp j P{u)du^ sel. (58) 

We integrate from a to t and obtain, 

exp J (3{s)ds^ z{t) < J (3{s)a{s)exp J (3{u)du^ ds tel. (59) 
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From assumption (l53l) and fl59l) we now derive the desired inequality, 

u{t) < a{t) + z{t) < a{t) + exp( [ (3{r)dr] [ (3{s)a{s)exp (- I (3{u)du ) ds (60) 



a{t) + J (3{s)a{s)exp (^J^ (3{u)du^ ds, t e I. (61) 



4.1.2 Bounds on the perturbations 5X 

From (H9l) we derive the following bound on the formal solution 6X 



/oo poo 
\Grit,t')5T{t')X{t')\dt' + / \Gr{t,t')5T{t')5X{t')\dt'. (62) 
-oo J — oo 

We will now use the fact that, by virtue of (152|) . where nonzero, Gr(t,t') = C{t',t). Hence 
(cf. ( !24ll ). inside the integral, 

\G{t,t')\ < ^ (I/ImIMOI + \fmh\M) = 9{t'), (63) 

with I/I A/ and \h\M being the absolute value maxima of these functions on the integration 
domain. The integration regions are in fact finite, since fl52p vanishes unless to < t' < t: 

\5X{t)\< t \g{t')6T{t')X{t')\dt' + t \g{t')6T{t')6X{t')\dt'. (64) 

•'to •'to 

Since g{t') is independent of t we can now apply Gronwall's inequality to obtain 



\SX{t)\ < / \g{t')6T{t')X{t')\dt' 

J to 

+ j'^ (^j' \g{t")ST{t")X{t")\dt''^ \g{t')ST{t')\exp (^j^' \g{t")ST{t")\dt"^ dt' . (65) 

On the interval (to)^) we assume the existence of a maximum of |X| and of |5T| and we 
call these \X\m and |5T|a// respectively. We also assume the integral \g{t')\dt' can be 
bounded by a number M. If 

-t 

\g{t')\dt' < M, (66) 



to 

then it follows that also 

ft 



/ \g{t")\dt" <M. (67) 
Jt' 



We thus find 

\5X{t)\ < \X\m {M\6T\m + M^TlLexp (M|5T|m)) • (68) 

The second term on the right-hand side is negligible compared to the first one for sufficiently 
small |5T|. 
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4.2 Solutions away from the singularity 

In regions I and III we will take 

T = + k/t\ <5T = i O (^^ j , (69) 

with h defined in The solutions to the unperturbed differential equation fHSl) are given 
by ^ 

^/\t\Jc{\nt\), y/\t\J^^{\nt\), a = a--, (70) 
where the Bessel functions, Jq,(x) and J_a(x), satisfy the differential equation 

5' . . d 



X 



2 



Mx) + x—Ja + (x^ - a^) Jaix) = 0. (71) 



dx'^ dx 



(This Bessel-negative-order-Bessel basis is more convenient for our purposes than the often- 
used Bessel-Neumann basis, as it approaches and for small values of t without 
mixing the two powers.) 

The unperturbed "compression factor" in region I is then 

^ u ^^ /rn /r^ Ja{-nti)J-a{-nt) - Ja{-nt)J-a{-nti) 
VlilVM w[^\U-nt),^\J^^i-nt)] 

Using the series expansion of the Bessel function for small arguments (they will be evaluated 
at t = — te), 

^°W~(i)°r(^- "^-1.-2.-3..... (73) 



we can estimate the various contributions to (1681) . thereby constraining the correction to 
the unperturbed "compression factor" . One can distinguish three cases: 

1) a > 1, Jo(— nti) 7^ 0, which yields 

|C(ti,g| oce(^-^)(^-"), \C\m oc e^^-'^^^-^l \6T\m oc e'"-^ M oc e^^-^^^^""). (74) 
From (J68l) . 6C{ti,t^) is negligible compared to C(ti,te) if 

c > -. 75 

a+b-1 ^ ^ 

2) a < 1, Ja{—nti) 7^ 0, which yields 

ICituQloce^'-'^^'-^l |C|AfOce°, \5T\m oc e^"-^ , M oc e°. (76) 

From (l68l) . 6C(ti,t^) is negligible compared to C(ti,te) if 

3 — a 

b+l-a ^ ^ 
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3) Ja{—nti) = 0, which yields 

|C(ti,t,)| oce(^-^)^ |CU./oce°, |5T|m oc e^^-^, M oc e°. (78) 
From (l68l) . 6C{ti,t^) is neghgible compared to C(ti,te) if 

c > 79 
b + a ^ ^ 

Whichever of the three cases is reahzed, it suffices for c to be greater than a number less 
than 1, in order for the corrections to the unperturbed "compression factor" to be negligible 
for small values of e. The discussion of interval III is completely parallel to what we have 
just presented. 

4.3 Solutions in the near-singular region 

The "unperturbed" equation in region II, 

X{t) + -Mt/e)X{t) = 0, (80) 



dt^ e 

is precisely that of the string center-of-mass motion. In order to simplify derivations, we 
shall assume a+ = a_, as required for well-defined zero-mode propagation (see section [3]). 
The unperturbed "compression factor" in region II takes the forno 

. , , Q22{\)\mY" - QiMU-^t) - Qi2(A)it.r-"4-v"'i + Q2i(A)it,rt^ei-2" 

CiAu.tj) ^--^ , 

(81) 

where the 2x2 matrix Q is defined by (1371) . and we have used iy[|t|'^, \t\^ "] = 2a — 1. To 
study the perturbation we will ffist perform the scaling transformation r] = t/e, Y{ri) = 
X{rie), which yields 

^Y{r]) + {e^r? + A^](r/)) Y[ri) = 0. (82) 

We now take 

T = Xn{r]), 5T = eV. (83) 

If we now choose f ^ rj"-, g ^ rj^^"- in fl63|) . M of fl66l) for the region {—t^j e.t^j e) (whose 
size, in 77, is proportional to e"'^) becomes (with the three factors coming from /, h and the 
size of the integration region): 

M oc = e-2c_ ^84) 



^This formula follows from ([29]), ([30]) and (l32|) via (l24|). It is also the same as (44) of pi. The expression 
given corresponds to small values of e. The corrections are suppressed by powers of e and do not contribute 
to the singular limit. 
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Because there are only power laws involved in flHTl) . the maximal value Cm is of the same 
order as \C{—t^,tf)\. Furthermore, |5T|a/ oc by construction. It then follows from fl68|) 
that 

|<5C,,| < (O (6^-2^) + O {e'''%xp (e^-^^))) \Cn\. (85) 

The correction is negligible for any c < 1. 

A subtlety in our above derivation deserves a comment (we would like to thank the 
JHEP referee for raising this point): one might have thought that the factor of in 5T 
of fl83l) competes with the smallness of e and undermines the validity of our considerations 
(for sufficiently large mode numbers). It is indeed true that, for each value of e (each fixed 
resolved space), our analysis is only valid for modes with sufficiently small mode numbers 
(though this range of validity increases infinitely as e is taken to 0). However, since the 
modes are completely independent, the limit for the motion of the entire string (if it exists) 
is exactly the same as if it were computed mode-by-mode. For that reason, n can be kept 
fixed in the derivations of this section, and the problem of competing with the smallness 
of e does not arise. (This attitude guarantees reproducing the e ^ limit correctly for the 
entire set of modes, though it does not allow to draw conclusions on the uniformity of this 
limit with respect to n.) 



4.4 Effective matching conditions 

Having analyzed the "compression factors" on subintervals I, II and HI, we can combine 
them into the total "compression factor" by applying fH3|) . As has been shown above, 
there exist a number c in fB51) between and 1, such that the "compression factors" on 
subintervals I, II and HI can be well approximated by the simplified expressions fl72]) and 
(IHTl) . with corrections suppressed by positive powers of e. One can then substitute (!72|) and 
( IHTl) into the right-hand-side of (l43l) . 

For a > 1 {Xk < 0), the Bessel functions featured in (1721) blow up near the origin (the 
inverted harmonic oscillator is propelled off to infinity). This threatens the existence of 
an e — i> limit. In appendix A, we display the divergences arising for a > 3/2. (For 
1 < a < 3/2, the limit may exist for individual string modes, but a consideration along the 
lines of section [5] would still indicate no consistent propagation for the entire string.) In 
any case, we shall not explore this case further since, as will be explained in section [6|, free 
strings are not a good approximation to motion in such plane waves. 

For a < 1 {\k > 0), substituting ([72]) and ([HI]) in (US]) yields 



C{ti,t2) = — [Q22WJa-l/2{-nti)Ji/2-a{nt2) - Qll ( A) Ji/2-a (-?^t 1 ) ^a-1/2 (^^2) 

+Q2l{>^)^^''^''Ja-l/2{-nti)Ja-l/2{nt2)'yn 

-Qi2(A)e2«-i Ji/2_a(-r2ti) Ji/2-a(nt2)7;:') , tl < 0, t2 > 0, 

(86) 

where 7^ are numbers originating from the coefficients of the power law expansion of the 
Bessel functions. 



14 



Note that the expression fl86l) has the same algebraic structure as the one derived for 
the center-of-mass motion in |9|, except that |t|" and are replaced by 



t| J_a(|t|). Demanding that the e — > limit should exist results in the condition 

Q2i(A) = 0. (87) 

It is exactly the same condition as the one for the existence of a singular limit of the center- 
of-mass motion (generically leading to a discrete spectrum for A). Under the assumption 
of fl87|) we obtain in the singular limit 

■ Q22WJa-l/2i-riti)Ji/2-aint2) - Qll ( A) Ji/2-a( -^^^1 ) (^^2) 



ti<0,t2>0. (88) 

The matching conditions across the singularity can now be derived rigorously by construct- 
ing two independent solutions to (139|) . Note that all the information necessary for such 
construction is encoded (cf. (HT!) ) in the "compression factor" given by (IHHI) . A convenient 
shortcut for this procedure is to recall the representation of C{ti,t2) in terms of two 
arbitrary independent solutions f{t) and h{t), and to read off the corresponding singular 
limit of the two solutions directly from (IHHjl . Writing Qii(A) = q and Q22W = —^/q, we 
obtain as a basis of solutions, 



n(t) = V=tJa-l/2(-^t), Y2{t) = y/^Jy2-a{-nt), t < 0, 

Yi{t)=qViJa-i/2{nt), r2(t) = -— Ji/2-a(nt), t > 0. (89) 



5 The singular limit for the entire string 



As we have seen in the previous section, for kX > 0, consistent propagation of the string 
center-of-mass across the singularity guarantees that all excited string modes also propagate 
in a consistent fashion. This is not sufficient, however, to define a consistent evolution for 
the whole string, since even small excitations of higher string modes can sum up to yield 
an infinite total energy [1]. As we shall see below, the condition of finite total string energy 
(after the singularity crossing) turns out to be very restrictive. 

The total string excitation energy can be conveniently expressed in terms of the Bo- 
goliubov coefficients for the higher string modes. To compute the latter, we shall form two 
different bases of solutions from flS^ corresponding to purely positive and negative fre- 
quencies at large negative and large positive times. More specifically, using the asymptotic 
expansion for the Bessel functions 



we construct 



J±a{x) 



sin(a7r) 




2 / TT 71' 
cos X =F I : 

nx V 2 4." 



X 



oo, 



— exp(iQ;7r/2 — in/A) exp{—iaiT/2 — in/A) 
exp(— '?a7r/2 + in/A) —exp{ian/2 + in/A) 



Yiit) 

y2{t) 



(90) 



(91) 
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such that, 



Tin 



-exp(mt), 02 {t) 



-exp(— mt), 



-oo. 



(92) 



Analogously, we introduce 



gsin(a7r) 



exp{—ia7r/2 + iir/A) q'^exp{ia7T/2 + in/A) 
— exp(iQ;7r/2 — Z7r/4) — g^exp(— ZQ;7r/2 — in/A) 



y2{t) 



(93) 



such that 



01 W 




Tcn 



exp(mt). 




Tcn 



exp(— int). 



-oo. 



(94) 



The two bases are related by a matrix made of Bogoliubov coefficients «.„ and 



[0^1 




an 


Pn 




'0r" 






P*n 


< 




02. 



(95) 



For the Bogoliubov coefficients, we obtain the following expressions, independent of n: 



Or. 



f3n 



2q sin(a7r) ' 
.exp(— ivro;) + g^exp(z7ra) 
2gsin(a7r) 



(96) 
(97) 



Here, a = \/l — AkX/2. The total mass of the string after crossing the singularity is given 

by [1] 

M = Y,n\l3n?. (98) 

n 

Since the /3„ are n-independent, M can only be finit(0 if /3„ = for all n. For k\ > 0, this 
cannot be achieved, since < a < 1/2 and q is real. (For A; = 0, which is the case of the 
"lightlike reflector plane" of [Hj, all /5„ will vanish if = 1, which is satisfied automatically 
for any reflection-symmetric f2(A).) 



6 Discussion 

Before we recapitulate our main results, it shall be appropriate to make two observations. 

^In general, one needs the uniformity of the e — > limit of /3„ with respect to n in order to analyze 
infinite sums as in (|98p . As remarked at the end of section our considerations allow to draw immediate 
conclusions on the existence of the limit, but not on its uniformity. However, since M is a sum of positive 
numbers, it is obvious that it will diverge when the /3„ approach an n-independent non-zero value (in 
the e ^ limit), irrespectively of whether this approach is uniform in n. For that reason, no further 
considerations are needed to draw our conclusions. 
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First, one can ask what kind of cosmological singularities gives rise, when the Penrose 
hmit is taken, to the plane wave singularities we have been considering. According to [7], 
if one starts with isotropic homogeneous cosmology of the type 

ds^ = -dt^ + e^dx'dx\ (99) 

and performs a Penrose limit, one obtains a plane wave of the form ([TH2]) with 

h 



k\ = , (100) 



Thus, positive values of kX correspond to positive h, i.e., Friedmann-like Big Bang singular- 
ities, and negative values of kX correspond to negative h, i.e., an infinite-expansion rather 
than an infinite-contraction singularity ("Big Rip"). 

Second, the dilaton field (discussed in more detail in appendix B) in the backgrounds 
of the type ([IP takes the form [2] 

= 0o + cX+ + ^lnX+ (101) 

If kX is negative, this expression blows up near (and so does the string coupling) posing 
a serious threat to the validity of perturbative string theory, and of free string propagation 
as zeroth order approximation thereto. 

For this reason of limited validity of the free string approximation when A; A < 0, we 
have paid relatively little attention to this case. What we could see is that, generically, it 
is hard to make excited string modes propagate consistently across the singularity (though 
it may still be possible to arrange such propagation by means of a judicious choice of 
the resolved profile ^(77) of the plane wave). The issue, however, cannot be competently 
addressed within perturbative string theory on account of string coupling blow-up. Our 
considerations can be seen as a motivation to study these backgrounds in the context of 
non-perturbative matrix theory descriptions of quantum gravity (the Matrix Big Bang case 
of |1] corresponds to the 11-dimensional analog of the plane waves we have been considering 
compactified on a light-like circle with k taken to —00 f5].). Some steps in this direction 
have been taken in |5j. (Alternatively, one could try to construct plane wave backgrounds 
of the type where the curvature of the metric is compensated by non-zero p-forms, 

rather than the dilaton, thus avoiding the dilaton blow-up problem.) 

For the case of positive kX, i.e., those plane waves that arise as Penrose limits of 
Friedmann-like cosmologies, it turns out that individual excited string modes propagate 
consistently across the singularity, whenever the center-of-mass of the string does. In those 
cases, the dilaton fllOlj) is actually very large and negative near the singularity, and one can 
expect that free strings are a good approximation as far as propagation across the singular- 
ity is concerned (the string coupling is small in the near-singular region). However, for free 
strings, we find it impossible to maintain a finite total string energy after the singularity 
crossing, provided that the (scale-invariant) singularity is resolved in a way that does not 
introduce new dimensionful parameters. The only way out appears to be to allow hidden 



17 



scales buried at the singular locu^ (even though the space-time away from the singularity 
is scale- invariant). To contemplate the possible physical origins of such dimensionful scales 
is an interesting pursuit, outside the scope of the present publication. 

Another relevant consideration would be the propagation of strings across plane wave 
singularities stronger than 1/(X+)^. Unfortunately, at present, little can be said about this 
case, even for the center-of-mass motion. 
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A Divergences for the case of the inverted harmonic oscillator 

As remarked in section for the case of kX < (inverted harmonic oscillator), divergences 
may arise in the evolution of excited string modes. These divergences may be seen via a 
blunt application of fj43|) . but it will be more instructive to make their algebraic structure 
more explicit. 

To this end, we shall derive a slightly different representation for the total "compression 
factor" in place of (H3l) . One can start by rewriting (!24|) as 

For any two sets of solutions {/, h} and {F, H}, the following relation holds: 

m\ _ 1 (W[f,H] -W[f,F]\ fFit)\ 
h{t)J W[F,H][w[h,H] -W[h,F]J [Hit)]- ^'^^^ 

One can then take four sets of solutions: one approximated by {\/^Ja{—nt), y/^J^a{—nt)} 
in region I, two approximated by {Yi_(t/e), Y2-(t/e)} and {Yi+(t/e), Y2+(t/e)} in region II, 
and one approximated by {^/iJa{nt), \/tJ-a{nt)} in region III. One can then start with 
fll02p written with the first of these four sets of solutions. In this representation, the 
functions featured in fll02p are easily evaluated at ti < 0, but not at t2 > 0. One then 



^If arbitrary resolutions, more general than (j4|), are allowed, for a given string mode, one should be 
able to reproduce (virtually) any matching conditions. This can be seen by assuming a particular form of 
solutions to the harmonic oscillator equation describing string propagation, and then reconstructing the 
plane wave profile necessary to produce this assumed motion. However, it is non-trivial to fit matching 
conditions for the entire tower of string modes in a particular geometrical resolution. For example, it is 
not obvious whether the matching conditions postulated in 2 should have any geometrical interpretation 
at all. 
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consecutively applies (11031) . fl32l) and (11031) again to insert the remaining three sets of so- 
lutions, with the Wronskians in (11031) being evaluated at the boundaries of sub-regions. 
In the resulting expression, all the four sets of solutions occur only with the values of the 
arguments for which we have convenient approximations to these solutions, and the total 
compression factor can be evaluated. As a matter of fact, this is simply another way to 
write (H3|). 

The divergent contributions to the total "compression factor" can be identified with 
particular Wronskians emerging from (11031) . when one constructs the total "compression 
factor" with the procedure outlined in the previous paragraph. For example, at the bound- 
ary of regions I and II, the following Wronskian occurs: 



W[V^tJ_^{-nt), F2-(t/e)] (104) 

The leading terms of both functions featured in the Wronskian are proportional to 
and therefore cancel by virtue of antisymmetry of the Wronskian. However, the sub-leading 
contributions have a different functional form and do not have to cancel. For example, for 
a > 3/2, one may consider the contribution from the first sub-leading power-law correction 
to the Bessel function, and the leading term in Y2-. This term is proportional to 

It]'-"] ~ (i05) 

and furthermore it is not accompanied by any powers of e in the total expression for C(ti, ^2)- 
For that reason, evaluating this term at t = — and taking the e — > limit will produce a 
divergence. 



B Background consistency and the singular limit for the dilaton 

As we have seen in the course of main exposition, consistent free string propagation turns 
out to impose extremely stringent constraints on the treatment of scale-invariant dilatonic 
plane wave backgrounds. For that reason, it was not crucial for our picture to explore fur- 
ther conditions arising from supergravity equations of motion imposed on the background. 
However, for methodological completeness, we shall present considerations for the singular 
limit of the dilaton field, and examine how this condition combines with propagation of 
individual string modes. These derivation will not have much bearing on the outcome of 
the analysis in the main text, but they may be useful for pursuing various modifications of 
our present set-up. 

If a time-dependent dilaton is used to support the curvature of the metric ([TH2]) in the 
context of string theory, the condition for conformal invariance of the world-sheet theory is 
given by [2] 

R^, = -2D^D,(P. (106) 

We shall impose this equation for all X+ in the resolved plane wave profile, and then examine 
the singular limit of the solutions for the dilaton. This is in contrast to the approach in 
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[2J, where the background consistency conditions at the singular locus were not discussed. 
The condition for conformal invariance (11061) leads to the equation 

m = (107) 

for the dilaton where d is the number of transverse dimensions X\ We want to consider 
the limit e ^ of the solution (p to this equation. In order for this limit to exist, the 
regularization Q will have to fulfill extra conditions. Since, in the singular limit, the space- 
time is regular away from = 0, we can construct a solution (j){t) to the left of the 
singularity and another solution 0(t) to the right. The requirements for the singular limit 
of to exist then reduce to demanding that the jumps in 0(t) and in its first derivative 
(j){t) are finite: 

h 

*2 Xd 

i<l>ii) + — I V^{v)dv (109) 

ti/e 



A(f)= I <p{t)dt= t0(t)J^ - / t(p{t)dt (108) 



il 



2 



and 



A(f)=- I ^Q{t/e)dt = ~— Q{r])dr] (110) 



2e^ " ' 2e j,^/. 



Thus, A0 can only be finite if 



-oo 

n{7])dr] = 0. (Ill) 

-oo 

If that is the case, the first term in fll09p is automatically finite, and we are left to demand 
finiteness of the second term 

lim / T]il{ri)dr] < oo. (112) 
If fl is even and satisfies ([5]), this second condition is automatically satisfied. 



B.l An explicit example 

We would now like to show that it is possible to combine the finite dilaton condition fillip 
with consistent propagation of individual string modes. Given the considerations in the 
main text, this translates into finding ^1{t]) such that fillip is satisfied and, in addition, 

—Yir^) + \nir^)Yir^) = (113) 

has a solution approaching ¥{7]) oc 7]^'°' for t] ±oo. We shall apply inverse reconstruction 
to fl{ri), assuming some shape of this solution and adjusting it so as to satisfy 

^ Y"(r,) 

-VrYdv = 0. (114) 

This "inverse reconstruction" technique is generally useful for contemplating qualitative 
properties of various plane wave profiles in relation to the singular limit. 
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B.1.1 No-go theorem for Y{r]) without zero crossings 

In constructing an appropriate Y{r]), it is important to decide whether it should have zeros. 
If Y has no zeros, Y'/Y is regular everywhere, and we can rewrite (11141) as: 



+00 



00 



Y{v) 



drj 



hoo , 



We now use Y{ri) oc rj^ for 77 — * ±00, yielding 



Therefore, if Y has no zeros, it is impossible to construct an ^(77) that integrates to zero. 
One must permit zeros (say Y{rii) = 0), and it is necessary to have bending points {Y"{rii) = 
0) at the same locations due to flll3p . We will aim at constructing a symmetric Q, assuming 
that Y is symmetric and restricting our analysis to rj > 0, and we will look for Y that has 
only one zero for 77 > 0. 

B.1.2 Piece-wise construction of solution 

We will now prove that it is possible to construct an Q that integrates to zero for a Y that 
has one zero-crossing. Q can be made arbitrarily smooth but for the simplicity of the proof 
we will allow Q to have discontinuities. The main idea is to split the contributions to the 
integral 

POO 

^iv)dv, (117) 







into two parts, separated by 77 = tjm- The part 

POO 

/ ^iv)dv, (118) 
will be chosen to be always positive. Then we prove that the contribution 

n{ri)dri, (119) 





can be made equal to any negative number while keeping the r] > r]M region intact. There- 
fore the total sum f lll7p can always be taken zero by adjusting the rj < tjm contribution. 
We rewrite equation (11131) as 

1 Y" 

fi = -- — , (120) 
and we take a piecewise Y{ri) (with a continuous first derivative), 

y(7,) = r-- -mi<v<VM (^21) 

\v\ > Vm- 
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Figure 1: Piece- wise construction of Y{i]). 



The function Y2 is fixed throughout our considerations, and we demand that it asymptotes 
to the subdominant solution for large rj: Y2 'q^~°' with 2a = 1 + \/l — 4A. As mentioned 
above, because of the denominator Y in. VL there needs to be a bending point for each 
crossing of the 77-axis. Y2 /Y2 is negative everywhere at 77 > rjM- The splicing point rjM is 
taken to be a minimum, and we demand that Yi{rjM) = Y2{rjM) = ^(%/)- We take the 
following ansatz: 

YM = iC- F(r^M)) (4- - + C. (122) 

A pictorial representation of our assumed solution is given on Fig. [TJ Due to the piecewise 
construction of Y it is clear that J Q{ri)dri consists of a separate Yi and Y2 contribution. 
The contribution of Y2 (i.e., -f^^ Y2 /Y2dri) will always be positive. It remains to be proven 
that Yi can contribute an arbitrarily negative value for fixed Y^tim) and rjM- With rjM > 
and A > 0, this is equivalent to asking that 

^^'"^ r dy (123) 

y + c-y{^m) 

can be set equal to an arbitrarily positive number. We know that Y{r]M) < C* < 0, since 
Yi should not cross the //-axis and rj = rjM is a minimum. First, if C = Y{riM), the integral 
above is 0. Then, for C 0^, with 6 = —C/{C — Y^tim)) > 0, we find in the limit of 

^ dy =. (124) 



/o y^-2y^-S 2V26 

For C — i> 0^ or (5 ^ this becomes arbitrarily large and positive. As a consequence ( 1119P 
can be made equal to any negative number (between and —00), and (11 111) can be satisfied 
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by appropriately adjusting Yi{ri). 
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